
visual_plot_t1 =
  data_4 %>%
  ggplot(aes(x=forcing_final)) +
  geom_histogram(bins = 60, fill = "darkgrey", color = "black") +
  xlim(-10, 10) +
  geom_vline(xintercept = 0, linetype = "dashed", color = "red") +
  theme_bw() +
  theme(panel.grid.minor = element_blank(), 
        panel.grid.major.x = element_blank(),
        axis.line.y.left = element_blank(),
        axis.line = element_line(colour = "black"),
        panel.border = element_blank()) +
  xlab("Integrated Score Forcing Variable") +
  ylab("") +
  theme(text=element_text(size=10))



visual_plot_t3 =
  data_3 %>%
  filter(abs(forcing_skd) < 20) %>%
  ggplot(aes(x=forcing_skd)) +
  geom_histogram(bins = 39, color = "black", fill = "lightgrey") +
  geom_vline(xintercept = 0, linetype = "dashed", color = "red") +
  theme_bw() +
  theme(panel.grid.minor = element_blank(), 
        panel.grid.major.x = element_blank(),
        axis.line.y.left = element_blank(),
        axis.line = element_line(colour = "black"),
        panel.border = element_blank()) +
  xlab("Basic Competence Examination, Forcing Variable") +
  ylab("") +
  theme(text=element_text(size=10))


ggsave(plot = visual_plot_t1, filename = "./_4_outputs/figures/figure_a3b.png", width = 5, height = 3)
ggsave(plot = visual_plot_t3, filename = "./_4_outputs/figures/figure_a3a.png", width = 5, height = 3)


